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A theorotical treatment has béen developed for the problem of an electrically 
small loop antenna loaded by an infinitely long, homogeneous, isotropic but hosuy 
ferxite rod, The loop which carries a constant current has. been idealized to be a 
delta-function generator, An effactive magnetic current (volts) is expreased 
explicitly in the form of an inverse Fourier integral, The contribution to the total” 
current {rum the aimple pole which can be associated with the surface wave is 
calied the tranamissicn current while the contribution from the branch cut giving 
rige ta the radiated field is, correspondingly, the radiation current. Alpo, the 
asymptotic behavior of the curront very near the delta-function source ie-investigated 
Two: values of electrica) radil of the rod are conaldered and for one of the cases the. - 
magnetic current ln plotted for a range of values of the permeability of the ferrite roc 
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ELECTRICALLY SMALL, LOOP ANTENNA LOADED BY A HOMOGENEOUS AND ISOTROPIC 
FERRITE CYLINDER ~- PART I 


By 


D. V. Girt 
Division of Engineering and Applied Physics 


Warvard University + Cambridge, Magscachusetts 


ABSTRACT i 
A theoretical treatment has been developed for the problem of an elec 
trically small loop antenna loaded by an infinitely long, homogeneous, 1A0- 


tropic but loasy ferrite red. The loop which carriea a constant current hae 


bean idealized to be a delta-function generator. An effective magnetic cur~ 
. 


rent (volts) is expressed explicitly in the form af an inverse Pour Ler inte- 
gral. The contribution to the total current from the simple pole which can 
be associated with the surface wave is called the transmission current while 
the contribution from the branch cut giving rise to the radlated field ia, 
correspondingly, the radiation current. Also, the asymptotic behavior of the 
current very near the delta-function source is investigated. Tro vaiues of 
electrical radii of the rod axe considered and for one of the cases the 
wmagdaecic current ie plotted for a r nge of valuen of the permesbility of 


the ferrite rod. 
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1. INTRODUCTION 

This report addresses itsulf tn the problem of farrita-cored Loop ancen-~ 
nas. Circular ioop antennas with permeable cores have been used extensively 
dn radio receivers. More recentiy, the radiative properties of loop antennas 
with spherical ferrite coves have been studied both theoretically and experi-~ 
mantally by several researchers. Loop antennas with cylindrical ferrite 
cores have not been used as transmitting elements, possibly bacause of a lack 
of sufficient theoretical and experimental informacion, 

By way of introduction, it is uaeful to consider an historical review of 
this class of antennas. In the early yeara of radio, receivers (540-1600 KHz) 
amployed a flat coil of wire, usually mounted on a Elat aurftace of the radlo 
cabinet, as the receiving element. Since the coil was air-corad, its perform- 
ance depended largely on the number of turns, coil] area and Q. With the de- 
mand for compact sets, it became increcaingly aifficult to place large~area 
coila far enough away from the chassis and get appreciable sensitivities. 

Out of this need for smaller sets evolvad ihe idea of uaing high permeability 
material for an antenna core and an early work reported on this subject is by 
Kihn, Harvey and O'Neil] (1940). Their experiments involved a core of finely 
divided iron preesed with a binder which soon provad to be uneconomical be-- 
cause of the large mass of material needed for a smail improvement. So, a 
large permaability material with a low loss was needed and found in ferrites. 

Since their use in broadcast raceivers, ferrite rod antennas have ra- 
ceived only exedutonet attention. As transeitting eluwants, they. hava beon 


atudied more racently. However, moat treatments [1]~-[4] have bean for spher- 


ideal ferrite corea; an exception is the work of Islam [5] which treata a cy- 


lindrical ferrite core driven dy a constant current carrying loop. The for- 


milation in (5] consists of finding the magnetic vector potential Kw A, ex- 


plicitly in an integral foma. Some numerical results are also preaented for 
the radfated field and radiation resiatance at low frequencies of the order 


of 300 KHz. In contrast with the work of Islam [5], the present report is a 


direct boundary-value approach to find the electromagnetic fields avarywhere. 
For this purpose an effactive magnetic current has been defined and evaluated. 
At least in principle, the othar quantities of interest can be derived from 
the magnetic current distribution 1f it is precisely known. 


IL. ELECTROMAGNETIC FIELDS OF A LOOP WITH FERRITE CORE 


Figure 1 shows an electrically small, filamentary loop antenna of radius 


a, loaded by an infinitely long, homogeneous, isotropic but losay farrite rod 


Sener aa 


of the same radius. The ferrite wedium is characterized by y = Hg (uy + iu). 


= Egfer + ten) and ky “wre. The medium sucrounding the rod and extending 


eee eile en, oA eee Ree om 


to infinity 1a free space, characterizad by w g and ky * wine g: The 


or & 
radiua of the loop ia much lees than the wavelength he in the ferrite mediun 
so that the loop current os ia in phase at all points ond essentially a con- 
etant. Thus, the only gource of stack tombenatic fields in this problem can 


be represented mathematically by $158 (6 - a)é(z). Furthermore, there is azi- 
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muthal symeetry so that the field quantities do not vary with respect to the 


$ coordinate. An harmonic time dependence factor exp(-iwt) 1a implicit in 
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all ffeld quantities. 
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At thie stage a discussion concerning the ralevant f1:ld componente is 


uote Pesta 


in order. Islam [5) states chat due to the symmetry af the problem, only the 


¢ component of the magnetic vector potential k axists and then procrede te 


yarn ag gsm 


a find Eyy Hy and 4 through Ay» satcing ail other: ideld cowponente aqual to 
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t zero. In evaluating field quantities in certa.a antenna problema, it ia cor- 


venient to uge the componenc of vector potential that {se parallel to the di-~ 


A, raction of che current in the antenna, wary, A, dn the case of the dipole 
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GEOMETRY OF THE PROBLEM. 


FIG. | 
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antenna and A, for a loop antenna without the core. Is both these cases the 
parallel component ia sufficient to solve the problem completely. The ap- 
proach geverally adopted ie to set up and solve an iitegral equation for the 
current on the antenna. With the current distribution know precisely, other 
quantities of interest can be derived. 

There are a few sim’iaritiag between the conducting cylindrical dipule 
antenna and the ferrite-rod antenna. The dipole autenna ia made up of a wire 
of high electrical conductivity and ie driven by a delta-function voltage 
generator wnile the ferrite rod antenna consists of a material of high perae~ 
ability (magnetic analog of alectrical conductivity)” and ie dziven by a 


delta-function current generator. Practical metals like coppar, aluminum and 


brass have high enough conductivities to justify an approximation of vanish- 


ing f1lelds inalde the materdal of the dipole antenna and, if more arccuracy ig 
required, cheories do exiat for imperfectly conducting dipole aatennas [6], 
(7). While the dipole avtaunat seotien hae been set up and solved with an in- 
tagral equation, che loop loaded by a ferrite rod is a boundary~value problem 
formulated in terme of d/fferential equations. However, on the basis of phys~ 
ieal mecherisne, the ferrite-zod antenna can well be compared with che dielec- 
tric rod antenna. [8]. In the ferrite material, the magnetic Aisoied cnet 
aligned in the direction of the magnetic field giving rive to an effective 
magnet ization 4 whereas tha electric dipoles get rearranged in the dielectric 
wedium to give rise to a polerization B, Thies analogy will be diwcussed in . 
further detetl at ao later stage. 


Returning to the queation of relevant field somponants, the loop caJries 


th fete 
- Pextaability can be called the magnetic analog of electrical conductivity 
wince conductivity and permittivity con be repressuted interchangeably in 


a watarial or medium with complex parameters. 
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an azimuthai current and excites the magnetic dipoles inside the forrite med- 


Since theae currents on the antenna ate in che ¢-direction, a component of 
magnetic vactor putenticl paraliel to the currents, Aye is suffictenc to.de= | } i 
rive all the non-vanishing field components. This is basically the reason 
why E. = 0, Ey = 0 and H, : 0, and Hs Wy and E, must be detarminad by solv~ i 


$ 
ing Maxwell's equations, appropriately written for various regions and with 


suitable boundary conditions, In fact, thia procedure dues not require a 


ium which can be viewed 46 microscopic current whirla as shown in Fig. 2. 24 
current diatribution to be dofinad on the antennas; however, a knowledge of an 
| 


equivalent current distzibution on the infinite rod could perhaps be very 


useful in peadicting the characteristics of a finite rod antenna. It is } 


+ 
mathematical’y inconveni:nt to work with tra V x M currenta depicted 1a 


ne 

ee 
@ 
re 
a 


f 3 Fig. 2; hence, an equivalent picture given in Fig. 3 ie used and the magnetic } 
: eurrent denaity i ide defined. It is a volumetric current with apecific D and . 
i z dependence, which can be integrated ovny the crous section of tha antenna i 
F to obtain an equivalent magnernic current 1¥(2) (volts). ‘This current can be ) 
i derived if the electromagnetic fielda inaide the ferrite rod are know. It } 


f 
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da now necagsary, therefore, to determine these fields which are solutions of 


Maxwell'u curl equations: 


vx de -~b | a 


vxie dad | a 7 (2) 


ane MMB och LEE nmin 2: . 


Eliminatiag A from (1) and (2) givea 
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“ $15 800 - a)é(z), B, = E. = 0 wad 9/39 = 9, and using au expan- 
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AN THE FERRITE MEGIUM BY 
THE CURRENT CARRYING LOOP, 


FIG. 2 MAGNENIC DIPOLES EXCITED 
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gion in cylindrica}, courdinateu, (3) reduces to 


2 2 
a 14 ( 25, OF } a | 
spe ei k® - 4+ He | BE (py2) = mtu (4) 


By solving (4) for E, and veing (1), 4, and H, ean be obtained in cterma of 


¢ 


KE, as followa: 
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In ordar to cotve (4), a Fourier transform pair le defined as follows: 


Uficigsty Aiea ina Daa atl pela a So eal ae Tcl A ple Nl 


eo > 
Eo.) * [ EsQoynpe’** ax (6) 
é 
Ey(ese) © Sk sy Ete, ott 7) 
The Fourier traneform of (4) is 
are ae el ee eee | a 
ar. . ao [ mE Ym 4] E(p,&) = ~Lwpt 6 (o - a) (8) 


p 


Continuity of the tangentisl electric field requires that E(p,£) be tantdnu- 
vusi at pa (first boundary condition). The discontinuity in the tangential 


magnetic fieid tr the true electric surface current density, or 


G2) a wD 
% z 


wage -1g0d) where the auperscripte refer to regions I and Il 


aa ghoun in Fig. 1. With (5b) this becumes 
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The Fourier transform of (9) ia 


ot 2 = t 
[se 2 azo!) : PASa . aes | bets ie 
dp u dp a ally. ‘sti oOo 
(second beundary condition) 


Now the solution for (8), which satisfles the homogeneous differential equa- 


tion and la single-valued, fulfills the above boundary conditions and 1s weld 


behaved at po = O and infinity, is 
1 | er , 
z‘ (40) * AJ, Bes ~ Ep for 0 <p <a (10) 


#2) ¢, Eb) « wi) ( J Ke - ce p } for a < 


a 
G 
A 

3 


(11) 


Let Yo Ke 7 e and 4” /« - ee; X¢ follows from the application of 


the boundary conditions (see Appendix I) that | 


7 ha 7 
[ CO =e 0 
J, (ya) + ay, J) (y, a) “Lut ni (y A) + au yYoH ey a)} duu uate 
L*'h. Mi1 0 ay 0 Or "O 
This autets equation can be eo lved for A and B to give 
A- oc (yp) (12) 
Bn Asal) 1672/06) : a . (15) 
vies D(E) is given by 
Pe reese (2) 7 
D{eE) Oly To (yyeny “Cy 9@) ~ Yous Iy6ranif Save, (44) 


The substitution of (12) aud (13) in (10) and (11) gives the cvausiorasl 
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flelds in the two regions. Thus, 
EO (oe) % duvet) cy aya, (y,0)/008) ; (15) 
B) (oye) = duuarts oan? (¥9¢)/0(6) (16) 


Te application of the Fourlar inve¢sion forsmuia Rives 


: ie (1) 
fwyal, «HH (y,8)J i0Y p) ; 
(e) + 0 - 0 hehe LER 4 
Ey (p,®) ar, x ~ By aes dt (1.7) 
12} , Loyaly a J yyy CD 0) ite 
Ey (p,z) * = EL sae ae ce’ oe (18) 


Once the praceding integrals vre wvaluated, the electromagnetic field ta com-. 


pletaly. detarnined 1£ ure ia made of (Sa) and (SE). 


Ill, REDUCTION OF FIELD EXPRESSION TO THE CASE OF SINGLE: TURN LOOP ea 


IN FREE SPACE 
When ae 1 and * 1, the problem is aquivalent a that of a loos an~ 
tenna: in free apace. Hence, with a a 1 in (18) che Fteld ahoutd reduce 
to that of a constant current-carrying loop antenna, ln order to Bobneye 
and . 


this reduction, tha quantitiva Wee fey and |! become 1, l, k 


Ups veopactively. With these changes (18) becomes 
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irs 


da in agreewant wlth the cusults obtained by King [10] and Wait [31]. 
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IV, MAGNETIC CURRENT ON THE FERRITE ROD 
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i From eke knowledge of che diectiie field inside the ferrite rod, the } 
: wagnatic current fin tho antenuc can be tound. As pointed out in meceeen Il, 
© a knowledge of thio current could be useful in erder to predict the charac~ 
i tertatica of a finite rod antenna, The fortowlnr procedure ts adopted in i 
E | (1) Q) can be found eine 


é ia known, HE 


(5b); 2) piace che ferrite-medium is assumed to he homogeneous and isotropic, 


; , 
finding the current NAS 1) -since E 


ui”) (6,2) “ Ch ~ pH (p, 4) is easily found from which 
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The use of this formula for MO) (2) dn (21) yields 


(1) 
eH (y,,a) a 
1(a) tal, ~ Darhy9[ Wi Gar craee atte J 3g 6¥ 0) 140 ao ) ar,| 


Let the variable be changed su that x * ¥yp3 then 
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T(z) ~-Lu(, ~ Vattug| eer. aa 7 } KJ (x) dx hac (22) 


To do the x integral, the following (dentity {fe umed: 
a 
a ' a fe 
mJ (x) mls (x) + J Ox) i [x (<) 
Roth aides can be integrated vith respect to x: 
! xJ 9 (x) dx = xJ, (x) 
Therefore, 
ayy AY) 
H KU y(x) dx = xd, (x) ; = ayy 3)(y,8) 


When this result ie substituted in (22), the following expression iw obtained: 
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Thus, the magnetic current t&) (volts) in the ferrite rod in expressed 


explicitly in an inverae Fourlar integral form. The investigatrion of singu- 
laritiea of the integrand and auwaericel evaluation of the integral Form the 


subject of Sections VI and VIL, respeccively. l 
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Pgh the use of tha modified Beswel functions: 
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With the use ci the asyaptotic axpansatlona, 
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& Since the current On the antenna. is an even function of z, one can write € in 
H : i 
f the foregoing expression as |£| and take its cosine inverze transform to get a 
i * 
i 1 (z) aez-+Q. Thue, 
* 
T(z) ge oe Shs 5: 
we Lbaul uy ne gn |[z| + Zinite integrals (24) 
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‘This equation states that the magnetic current: has a logarithmic eingularity 


at the source and is similar to the isolated dipole antenna in free space ob- 
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tained by Wu and King [12]. a 
VL. TRANSMISSION AND RADIATION CURRENTS ON THE ANTENNA 


To evaluate the following integral, tae singularities of the integrand 


must be investigaced: vA | i 
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The total magnetic current I(x) can he thought of ag uw gua of 2 trang- 


rotsaco MES, 
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* * 
miesion current 1,.(2) and a radiation cutrent T(z). The contribution from a 


‘ b simple pole gives rise to the tranamission current and is uavsoclated with a ‘ 
\ : surface wave on the antenna whereas the contribution from the branch cut is i 
Mg i correspondingly che radiation current that maintaing the electromagnetic ; 
x fielde at distant points. 4 
a : Note that € = tk, are not branch points aincea che integrand remaina un~ i 
*) oe 3 
ut if changed upou adding 1 to the argument of Yh: Thus, € = ko are the only a 
' i branch polutsa. The poles of the integrand can bu-dntoanined hy aolving 

: DCE) = 0, which will be discussed in derail with reference cto Fig. 4. This 

e figure shows the path of integration, the pele docation and the branch cuts 


in the complex € plane, At this stage, for illustrative purposes, YY and oe 


are assumed real. 


cl <>, witha = Jeo - ken ~iy, and 


a) On the real axis, for ky < | 
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Ru VJ 0? ~ 4% miemdy s DCE) = daal (toa yn? (apa) ~ bp dhal Ciuay ts?) (apa). j 
0 0 0 i rd 0 
Introducing modified Beusel functions, D(&) ~ 0 reautres that faal)(oa)K, (fa) «4 
+ H. B01; (aa)Ky (Ha) be equal to zero. Since for real and positive values of b 
iy 

a and f, the modified Bessel functions Io le Ko and K, are all real and i 


positive, this requirement cannot be met and hence uo pole can exiet on this 


part of the real axis. 
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b) On the part of the real axis where 0 < [5] < k, and on the entire 
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imaginary axis, D(E) = 0 requires that 
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FIG.4 COMPLEX €-PLANE SHOWING THE SINGULARITIES AND 
THE PATH OF. INTEGRATION 
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The laft-hand side of the above equation is alwaya real for the range of € 
values being considered whareas for the righthand aide to be real 
t] e s ! 1 
J yg) ¥g (7g) V9 S¥ g@ Fo lr 9") should be equal to sero. But this ie a 
Wronskian and cannot be equal co zero. Therefore, there La no pole on the 


part of the real axis for 0 < je] « k, or om the entire imaginary axis, 
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c) On the real axis, for ky < lel < k,, the equation D(f) = 0 becomes 


"1" 


y;a ody (7,a)HO? (482) - 4 Aad orang (4a) = 0 
or (26) 


aJ (ya) K, (Ba) + u.8aJ, Cy #) Ky (Ba) = 0 


“WJ gly) /5, OG) . HW aK Ox) /K, Oo) . (27) 


where x and y are both positive and real with x @ fa gad y = ¥yas The trans~ 
cendental aquation (27) de similar to the one obtained by Sommerfeld [18] in 
the problem of wavea on wires. However, the graphical method used here for 
solving the aquation io ensentially the same ‘as that of Duncan [8], since 


the riesbt-hand side of (27) is always positive and real for luealess ferrite 


medium, a solution 1a voasible only when y ia auch that Joly) and Jy (y) carry 


opposita signs, This can alsa be chsarved. tn Fig. 5 and leads to 
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whore £, 18 the solution, d.a,, D(tE,) = 0. 


fa ae | | 
<a/k tr ey. for i+ 1, 2, .0) |” (28) - 


a 
= 
3 
# 
= 


bed at a te ot! 


Tented FR Ler i igre ee Tio LOL dow ict OIL tego 


Se eel i “he 


Site Ce TRE 


Scr rc rr a ee ee ae ne ey : = “0m. 7 j 
Le et nn ea! ae ss 7 aoe ne ae Ness. Free IRB or eee ae 


‘ aes tae 
Pre 
4 


GEWDERTAL 0M (26) 


S 
M4 


‘ ort d 
t = —— o. “ee 


abe 


CAP / (Ap ao ey 


| 


GRAPHICAL SOLUTION FOR TRAK 


fae eer Me 
Ws 


Fig § 


20|--—-— 


"ois ble a ATM SES we, Sin SUGIMOTO tn ens Teal oc mae, eae nee oe oo 


we yeh : 
You 1 aaro of IgM 


we (se1yeh P 
Yy (+1) (£41) zero of &) ‘ 


for exawple, 


¥o.y 7 2-405 Yy 9 7 3-835 
Yo2 7 5+520 ¥y.3.7 7-025 
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From Fig. 5 1t de aleo clear that for every value of x there are infinite 


values of y which satisfy the transcendantal equation (27). Rach of these 
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i 2 selutions corresponde to 4 rotationally symmecrical TE propagating mode on 
f the antenna. Fig. 6 illustrates the multi-valued nature of y, arising out of 
e the infinite branches of the left~hand side of equation (27). Each peint 4 
(x,y) on the dashed curves in Fig. 6 leads to a possibie solution tEy: Also, ( 
ef xe af et we and y “4 Je -e which leads to 
x? + y? = RB? (29) ! 


uhere RB = (akg) (ec, aie 


Since x ard y have to satiafy equationa (27) and (2%) aitmultcaneausly, 


there ave now a finite nusber of solutions ac exemplified by the cirele C,. 
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‘Tf R tea auch that 2.405 < @ < 5,520, only the dominant TE wode is supported 
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by the farrite rod. If R < 2.405 like on Gy» the antenna is below cut-off 
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for wil the propagating surface modea. Furtheruore, for practical ferrites 


aincy bat, 7> La B Wak, thua, one can reach a conclusion that ak, has. to 
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be at least 2.405 for the aurface waves to appear and additional modes are 


“supported 3£ MY de iucreased sufficiently. Alac, when a surface wave ia 
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present, ita propagation constant fy will lie between the wave numberu ky and 
ky of free space and the farcite medium, respectively, 

d) it gtdil vemadne to axamine the finite complex € plane for che golu- 
tions of the equation DCG) = 0. Since no analytic solution in pussible, au- 
merical procadurea consisted of computing thn magnitude of the recivrocal of 
D(E) at sevaral grid potiuts in a square of size ky dn the first quadrant of 
the complax — plane, This computation was cartied out on both the Redirann 
aheety of the complex incegrand of (23) for Ue 10.0 and Bk = 0.0b. A 
sulution to D4) © O da taenndified with a peaked behavior of Increasing am- 
plitude f4 lL/p(é)|. ‘Lhe veal axda solucion a7. of ¢) with ko “ | 1 ok 
was foun! wiehin computational accuracy and ne other solution could be found. 
Although thie da not a conclusive search for the voota ef D¢f) = 0, db should 
ba Yermbas bub that the goluttons, Lf any, dn the lover half plane, leading 
to growi:e wavea, ara co be discarded. Furtherwore, any golutdon away trom. 
che taal arta in the upper half of the complex © plane gives rige to rapidly 
attenuation surface wives which are significant only a¢ vary short d4utancaa 
fyeom the delta~funcsdon BaNeLALOT . 

Keturning to the dutegral of equarion (23), the original pach of fute- 
gration C, whieh rune along the real axto with guitable indwatations, can be 
deformad and ehown equivalent to the contour I if the pole contrdbution at 
fom bo fa aultably caker fnto account (ace openken Tr). 
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The total magnetic current T(x) can now be written as 
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: : Tha contour [ as shown in Fig. 4 conuista of paths Gy and on which xun 
te he left and right of the branch cut and alao GC. which fa a se@ad-circular 
: path cround the branch tip. It can be shown thac the contribution at the cip 
c ; 
ee he da vanishingly small which leaves only sections Oy and Cc. to be computed, 
i . 
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Substituting far by apa oy from (32) into the aboyea and uring toe ana- 


lytic contimuakion propestige of Beamal and Hankel fudctdons, one obtalus 
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Es Upon using (32) and (34), the cwo integrals dn (33) can be cembined to : 
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wheres to,” 1207 ohme te the free apace characteristic impedance; 4(u,v), 4 
Biu,v) and u,v are datined ta (35) and (36), respectively. 
VII. NUMERICAL COMPYTAT*.ON 
A. Raddation Fart of Magnetic Current 
In oxder to compute numartially re 4e) from (37), dt in uaeful to examina 
the nature of the dutegrand. Equatlou (37) te reweltten de follows: 
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for coaples ‘anusaeaty iu uedge by ohicdn ify). Hubweue dma : 


BSLSML (of Bhet [15]) haw beer modified for double precdaion accuracy. This 


ie progxvan uses the series axpansion for JQ) 
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i wined by sulving a quadratic equation. Then uge is made of Horner's algor— ‘ 
f ! 
‘ ithe by caating the .»Lidte sertes in the form (2(z(2 + a) + b) + co)... and : 
cha productw exe computed frow the innermost to outernogt, thus idndm zing 
4 the round -off errora, Vor large atyumenta (|x| > 20), His: aaenranks expan 
: pion le uaad. — i 
Subroutine BEGH cvmputes thea Hankel funce lens ins ey Te aut We : 
; for complex argument x with double prectaiua accuracy, It maces uae of the : 
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| oe _Pragy was 'RESY.' rom the 08/360 16M Sedentifie Subroutine package which was 
modi fted by Dhar [ua] for complex argumenie. *BESK.' computer the modified. 
i f e Wenaal fuantton Ko Ce) and KO) which ave then uged ip computing the Hankel 
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paris of the total integrand are rapidly decaying functions of y and Lt waa 


found tat f. and re need be calculated for y ranging from 0 to 50 only, 
b , In the expression (39) for Ely) u and v are the complex crpumeuta of 
I; is. Bessul and Hankel functions, respectively. With: the maximum value of y noar 
{ By 50 and for tha range of values of ak,, »_ and ce, considered, ful and |v; ‘a 
not exceed 15 and 10, respectively. This euaures that the subroutines BSLSML 


and BESH are used wel! within the range of their validicy. 


1) €(y) as a. function of y, and the numerical integration for ti radia- 
tion current: 


In this section the behavios of the real and imaginary parta f. and t 


of f{y) (which can now be calculated using BESH and BSLSMLJ ta Tinucawens " 
fhe dieslactric coumtant fe of the ferrite madium te held constant at 10.0, 
Twn veiuea of olectrical radii, viz., aky = 0.05 and 0.1, are conaidered. ; 
Foe each aky the relative permeability We dn varied over u range of values | 
axceuding from 10 1 200. f. and fy are shown graphivally in Figs. 7 and 6. 
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For the rod with the spaller cadlua, t. and fi hava, reapectively, a iy 


positive and negative peak (Fig. 7) initdally, but aa Wie is increased their 


rOlee ara reversed. A gomewhat simntlur behavior iw found for the Larger 


radius (Flg. 4). Furthermore, in either cave, both i. and fy tend aaymptoti- 


cally to wero for large values of y, The decay of both the real und imag 


inary parts of the cotal dtutogrand (f(y) oxp(~yk,z) ] ia even faster because of 


the wultipiicative veal exponential factor. Due to this, a preliminary eval- 


uation of the integral of equation (58) showed that che upper Limit of inta- 


gration can be replaced by 20 or leas without any st{guificant lows in accur- 


acy. Lt ie important to perform the incegration accurately arouad tha peak 


because of ite significant contedbution te the total dutegral. A 12-point 


Gausn quadrature routine from 08/360 IBM Scientific Subroutine package has 
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bean used [16]. In order to meet a specified convargence criterion, tha 
total cange of intagration was divided into gsuffictent number of panels, noe 
to exceed 5 in any case. Since tha integrand has been previously calculated 
and plotted (Figs. 7 and 8), the iecation of the peaka {mn both tha real and 
{fwapinary parts ate accurately known. Panels are of unaqual width and are 
worn <\osely spaced around the peak. The optimum number of panels M ln de~ 
cided by requiring that the value of the integral using M and (M # 1) divi- 
sions differ by lean than 10" 1u wagnitude. The resuits of thease computa- 
tions are shown graphically in Figs. 10and 11; they are discussed in Section 
VIL. 
B. Transmiasion Part of Magnetic Cucrent 

The transmiedion current on the antenna td ativan by the contribution of 
the residue at the pola &— « bor to the integral cf aquation (25). This was 


calculated in (30) to be 
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The location of the pole § # £0 was briefly discussed in Section VI. Lt 
de now us#ful to sat up 4 graphical procedure to deteraine fy ard carry out a 


sunple calculation for the case of real parameters yu and fo. Fig. 9a) 


r 
slows the elecrrical radius’ ak, an a function of frequency renging from 1 to 
1000 MHz. Practical values of the diameter of the ferrite rod ars conuidered 
and it ranges from 1/2" to 4", In an actual experimental setup, care muse be 
taken to ensure the validity of the constant curvenrl eppvoximation in the 
dviving loop by requiring ako < O.1. Having datarsiued aky aud knowing ae 


and t.» one can obtain the va.ua of the parametes R wiich then ds plotted du 


Fig. 9%c) a8 tllustrated. A knowledge of Xy from Fig. 9{c) iu uped in ¥(d) 
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Co determina the value of the propagation constant bo normalized to the free 


tha dencainator and after simplifying, (40) bacomes 


| 

| 

{ 

2 space wave nuwber kgs Letting Gg/ky =n, carrying ove the differentiation in 
z | 

{ 

| 

| 


, ink,< 
‘ ; 
1. (x) “Angas” 

- om . (volta/amp) nee: Oo ae eet een coer ere ren cern caer ene s 

: fH 1 ¢ 

fs to eee FON eos “ooo 

a : 2 2 2 4 
; e ux, Fo (y,) xu, = 1) I ¢y,,) ri 
: og ro -1°°0 o'r 1-70 (41) 

2 where ug 120” ohme is the free space cheractwurdutic impedance and R * : 

: 20 21/2 zi 

y Tha traveling-wave nature of the tranemisalon curcgent can now be sean : 


from (41) a0 that it ia sufficient to plot the magnitude of the normalized 


a 
currence [1,.(2)/15] aga function of the normalized dietance. 


Exawupiei Wy 100.0, 6 10.0, Fm 50 Miz, 2a + 0 0762 u ov 3 in, 
4) From tig. 9(4), uk, = 0,04 i 

j 

11) From Fig. 9(b), K = 4.0 j 


A141) From Fig. “(c), = 2.50, Yo 4,12; thas, Xy/aky = 62.5 


0 
dv) From Pig. 9(d), nhs bo/ ko » 62.503 

Using the above values in (41) | 22) /78| tea found to be 39.05 elealans and 

ia show plotted in Fig, 1]. 

VILI. SUMMARY 


An electricaily exall loop that cacrles a conatant current end iz loaded 


eR ie a 


by an infinitely long, fou-genaouz, isotropia ferrice rod cas heen called the 
farrite-‘od autenna. The ferrite-red antenna is treated using a boundary- 

value epproach. An explicit axpression fer the magnetic current ds she form 
ot an iuveres Fourier integral laws bean devived and nuaarically coupuced, aa Pa. 


Two values of the aluctrical vadiue tux the loey ace considered. For wie of the 


cases the aaqguetic cusrent ie represented grephically ag a function of the 
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oraalized distance fcr a range ot valuga of the reiative permeability of the 
ferrite rod. The magnetic current ig found to consist of e@ trangeission and 
l He & radietion part. If y. and ec. of the ferrite rod are assumed to be reai, 
then the traummisezon current can be associated with an unattenuated TE sur- 
face wave. Thin surfaca wave is rotattonally symmetrical and has a cut-off 
condition. Since this surface wave does not contribute to radiatioa, the 
" ro cut-off condition 1s easily met at frequenctes where ferrite-rod antennas are 
useful in practice and, thus, the propagating surface mode can be made to 

a disappear. The redlation current on the ferrite rod 1s a decaying function 


of dfatance away from the delta-functton sourcs, Furthermore, the asymptotic 


o behavior of tha magnetic current near the delta~function generator was found 
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to be loparivhmic and, hence, similar to the elerttic current in the dipole 


antenna (Ws and King (22]). The analogy between the ferrite-rod antenra and 


the wouducting cylindrical dipole antenna was discussed in Seerion IL, It 


Was &lso moentioand that a comparison cf the ferrite-rod antenna with the di- 
elaciric rod antenua is possible on the basia of physical mechanisms inside 
the meacerfal., The preseat foraslacion can be compared diractly with the work 
vf Ting [17) on the dlelectric-soated aucenna. In this a currant diatribution 
which also includes « transmisuion and a eevaudine pate has baen obtatned. 
The magnetic currents plotced in thia report are for farcite cores of 
Anfinite leugth. However, in practice, low frequency anteanas Like the fer- 
rite vod are of necesaity fluite and even electrically iste: Therefore, 
logical extenadon of thiga formulerion iv to obtain wagnecic current distribu- 
tioia on @ Einite rod. With thiw current distribution pruciasly know, in 
principle, other quantities of iurerest like the radfsted field and radfacton 
efficiency can be derivad from it. If ie capected thac thie will form the 


aubject of Part {1 of thie report to be published at a latar date, 
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(Above functiona when primed mean derivatives with respect to their arsemants) 
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LIST OF SYMSOLY 
Gircular cylindrical coordinates 
= ig (us 4 4) parmeaahilityv of farrite madium 
Compiex vealative peurweabiiity 
- Fler + dey), parmitcivity of farrite medium 
Cowplex dielactcic constant 
Conductivity 
Strength of constant curwent in the driving Loop 
Radiue of loop ~ radius of ferrite rod 


Anguiar frequency 


z-transformsed siectric field } 
est : 
Seasgel function of firnt kind and order n 
Reuwsaun functlon of urder nu } 
Rankel function of firet kind and order n | 


Hankel function of second kind and order u 


~ w/e, wave number i 
Fourier tvanaform variable for z courdinato 
q 


7 
Yo ofy 7 | 
a eee} 
Vy ee: |. 
{(R, 4, 8) Spherical polar coordinates ' 
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DOUNDARY GONUTTLIONE aND EVALUATION OF CONSTANTS It YDELD BXCRESS LONG 


Tle purpose of thle aypendix 49 to avaluate the vonstants A and A owhhe 


appawsnud da equationa (10) aud (1), 
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APPENDIX II 
S <- THE CONTOUR OF INTEGRATION IN EQUATION (25) 1 \ 
ee: : 
: : Tha pucpese of this appendix da to sluplify the path of integration ap- ! : 
; pearing dn (25): : 
w© at) Cy adr, cya) ! 
: ® a 1 Li igieasae Léz ,. 2 
| i H (a) wodumugtu = DEQ fmm mga eo de : 
qi ue In the above equation the path of integration is the entire real axiu and is - 
Mi u 
/ te ; 
i. called the guouteur C, represented by the path A to N dn Fig. 12, i 
» Gonpsidaring the two closed patha AC)IUIPA and JKLMNOJ, 
i ! 
. OF if (+f oye fore (Li-1) ; 
ri A to it HP PA 
a: 
® } ()+ if (d+ ! () = ant (residue ac the pele nha) (11-2) 
: a to N yo OJ 
i leas y and tas ) are both equal te nero since the integrand is vanishingly 
& wall on the buge circle, Using tide resuls t- (1¥-1) and (11-2), 
ra 
é ! C)+] ©) Pe ‘a ~ f ( ) + Awd Crasidue at the pele me? 
; Ato H J to N uP oJ 
¥ 
8 
i Adding the gerd-ciraular path Wid co bor sides of thea above equarion, one 
a obkaduat 
a 7 | 
a [C) fC) + 294 (astdue ae che poke e=65) 
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This Ley the result used do squations (50) through (33). 
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APPENDIX Ii. 

Thim appendix essentially lisce all che Fortran IV programa that were 
ugead in the various computations. The ‘Main Program’ which appeara at the 
beginnadug was written te compute the radiation part of tha mapnatic current 
on che ferrite-rod antenna. Basically it involves a numerical {ntegration 
of a complex function. For thie purpose the behavior of tha real and imag- 
inary parts of the integrand for various parameter reacges was examined and a 
suitable Gaussian Quadrature routine was employed. The numericai evaluation 
of the integrand itself io comprised of cylindrical functions of complex ar- 


quaants. Previously avatlable prograns [13] were modifinad to maet the prea- 


ant. requirarenta, 
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c MAL PRUGRAM 
CHES RUST Ramey PAUECRALE COMPUTER THe AAUIATIUA PART OF TH MAGMETIC 
GC CUNARENE FESTA MUTE LAE Dt) AM PL EC TRICALLY SMALL LOOP ANTEMMA LILADED oF 
CA WIMGOENCOUS AMO TROT OPIC FRA LLTR CYR imei OF EMP IMT LENG The 
€ 
CU PUML TION BURPRAGAANE- 8 TRRLGTOL 
CC MPROUT INES § CAUSES, OF UM, OE Ney BSL SME ASV Z, ASTD ALVOQASVD ASV Og AMYE * 
& 
ooul CPL ICTY COMPLET ALT), REMLOR(TI 
acer CXYEMHAL TPL, TAS 
oes Comoe TUA, TEA, TGA, TRO 
O84 z RRAD (5,5) TUR, TER TROA 
G09 " FORMAT (21 DLOL40I Kd? 
aove AF CYUML@Y.9eN.n) GO TO 209 
aco? AGAD (9,20) TL. TU 
uove lo Fimmat(ztolasa,inds 
0009 Hee 
aaa WHIYE ¢e,20) 
vein 20 FORMATO (JORG *RAUTAT ION PAAT UF THE HAGMETIZATIOM CURRENT ON 
27KE FEANITE RUD ANTENNA F677) 
cor WRIYE (6,70) TUR, TERS TEOMA 
os FO FORMAT (90'S SK, 8 RBLATIVE PERALABILETY = 5 910.%,20R—f RELATIVE 
APURMETIVITY © FeDLD eae T TR RLECTREICAL AAUIUB KOA = "yULOee! 
ans Theti 
Gos Yaety 
aoe TheDFLOATIND 
\ GoiT WHITH (94) : 
" oole BO FORMATE AND LATS TH G MOK RAMS *RES! TL LOM eT IMAGT ROMO TAB? PLAN 
f CHPHRGEY LS) 
4 ao1y TUOL=TKOA 
! eos0 42 Yoe0.bd 
: oon, 40.06 
A ooar LGASTROL/THOA 
j 00?) TUPFLedya~ThISTN 
‘ oo1s 47 COMTI UE 
j vors CALL GAUBAI TLE TUSTHAR, THEA) 
I ude CALL BAUSECTUeTU, TAL IMAG) 
: aor? YyeTsoTREA 
i 9028 Tout beT 1 had 
4 ou29 Tye 
} 9030 TUSVUOCTasTLIZTM 
: ound TEP TUSLEAT2ZE WO YO 48 
0032 Thm (hSar-F23/2~00 
‘ Otay Timer 
H outs The ti t7o 
‘ oo3s 0 CONC TNUE 
' Upye GALL GAUSTITL STS TEARS TRIAD 
} ooar CALL GPUCAT TL TUTE as FIMAGD 
una TOY DOTAES 
! ons Veet aer (van 
04g Tisty 
9041 Tu=Tuet0 
{ Hoe2 1P (TULA .AS.0C} GU TO +O 
' does GLP MPLELTI Toh 
\ 0044 Fhe 19480 O01 0. 00, Led OCT MP1 Ue Oe) OOROTKEL) 
vas TADE@LOABRI TT) 
i eoee TREATREALIDE) 
| ons TIMAGSAINACLOID 
| posa TPH ELE, V7 De AGDSVZOT OOD FOATANT( TIMADSIRGAD 
‘ o0a9 VAITH C640) AGAQTNOL, TALA TI maa, TABH, CH 
£ ocxa 92 FORMAT CAS PRR AK DAG ear hy IML e MED?) 
? Ona) Thera 
t nose TUT 
: onsy TAQL©YNOL eT FOR 
} 9034 JP PORE N91) BU TO 42 
‘ Oo%S eo 2 
H ous 200) CONT Emut 
i oosT thy 
t 
t fC FAH OO4EOFOS 04 8 SOE BE SE Se EASES PE SETEE JES ae BEST EURO REEL YE Oee TEE 
H € 
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: MES PLOGRAM CITRPUDES THE INYEGRAL OF FCYRSENPI~YRg dd OCTWEEN Ti) 
; c SPECEFICU LIMSTS VIEMG & LE-POTHT BALES QUAUMATURE FLARES AS 
H C SINCE PLY) 18 & COMMER PUCTIUN, THE Stak amt [MAGINASY PORTS oKE 
: C LATECAAVED SEPARATELY. 
4 c 
i Cre y SOBRE IME PAUAST TE eTUSTLA THE ALD 
i aoa? TPLICIY KEALeR(T) 
4 ‘ Quus Tae POOEVIETLD 
‘ ea04 Tete ty 
, aoes TCH eh POTOOSE TIS DOSER WU DFTA 
4 9008 TABAL =. 29S CY COS OSTASPAAD= OU TAMA TROTCDATENLT ATC)? 
i eeor TU, GUZO9O07 31 847174 5008 18 
¢ oeds TRRAL WT AAAL © LS IACHSSZ ET ASEZ FU LOL Ti Md Fae tap oTand YA-TUd) 
t e007 Tew, SOePhE DS POU TINZ S40 U9 TH 
: serc TREALH TERE 2 OD ONTL ASR TIOTELAD- LOU PLR TAT YC UO TLMETA- TERS 
i 201) Ylon FO baBE+ FTL ARSC 209878 
¢ Goiz TRPALSTROAL OLOLSED TIT DO TIES EOSSITLALTASTC DOTA MITA TCDD 
i eos Thee hZ AVES F395 9G TEE CONS TE . 
: urs FRED CTREALO SP AGTOMEGBZAGA TIA EE (LRU TASTOD PTLME TART OD 
‘ ok Thm eB 2OUK7OSIS$ M345 BOL 8 Te 


eoia TAR AG=T Ret Tmbal ©.) 1937 35.2708 TGR OVOSE TLR TACEC beTLKETAYTC}I) 
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€ CDRS TRS y SESE 62E6 AE PERS VERS TE BST FE FOR SOS TE BEES OF FE OS be oT OSE O44G4 O49 : 
i ‘ CG OTHER FUMCTICN LUFBPROGRAM COMPUTES THE FAM. PART OF THE COMPLEX ci 
3 C RNVTEGAARO FEV) POR A GIVEN Y USING THE SUBAOUI EN *§ DFU F : 
tt Cc 
a : opal DOUBLE PHFCIATON PUbC TION TERACTYD 
iy a Qguz TMPLICEY amaeaeT) 
i os 00% CUMMOM YUM THR, TROAS TROL : 
A oe 40a CALL NPUML TUR THUD TEAS THOAS TY SYA! STEP 
4 cous (PRL a TRE : 
‘ Cyrey RETURN 4 
‘} 00y7 BRU 
Ve 1 . 
} Z Co FRRAHSL AES AETESE DESC ES ORES OE9 SSH ES! CREDEDRSS BESS HESS SE EEE DENTS SEE ERED : 
{ z a c i 
1 Cc i 
hI i C THIS FUNCTION SUBPROGRAM COMPUTES THE IMAGINARY PART OF THE CUM LEX 
4 i GC OTRTRGRAND FAY) bum A GIVEN Vo 
ree ¢ i 
4 ! & oyot COURLE PRECISION FUNCTION TREQTY) | 
t, an aooa IMPLICEY MEALOQ(Y) # 
a) ee aaoy COMMON TUR STCK, TROAYTROL : 
} a oou4 CALL MPUNETURS TROL TERS TROA, TY AVARSTIOD 
y : rer) ThLeT te i 
! £ 2006 RETURN 4 
" i @oo7 CAD 4 
hy : iw 
5 ® 4 
f u ‘| 
i CORNER OOO AMOR EO SERU THY RHO OER OF ORE OEE DORE N EET TF OF EY IT 88 FORESEES ‘ 
: ‘ c : 
: t ¢ 4 
y 6 BUBAOUTINE © OFUN PL 4 
C OTHIA SM ROUT He COMPUTES YH! COMPLEM IMTEGUAND FLY) FOR A OIWAH Ye j 
i C OTHE CUOMFUYED RESULT OF THIS BUMRCUTING 18 USEU Te THe FUNCTION 7 
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